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Abstract 


An experimental and analytical program has been carried 
out to evaluate the performance of sound suppression techni- 
ques in ducts that produce refraction effects due to axial 
velocity gradients. The analytical program employs a com- 
puter code based on the method of multiple scales to calcu- 
late the influence of axial variations due to slow changes 
in the cross-sectional area as well as transverse gradients 
due to the wall boundary layers. The overall objectives of 
the study were the verification of the analytical model 
through direct comparison of experimental and computational 
^results and the analytical- determination of the influence of 
axial gradients on "optimum"- 'liner properties. However, a 
weak non-uniformity in the standard parallel-duct . theory at 
the "optimum" liner conditions creates a singularity in the 
higher-order perturbation results (that is, the non-parallel 
effects) ; thus the analytical studies have been unable to 
examine the influence of nonrparallel ducts. on the optimum 

t 

liner conditions. However, /for liner properties not close 

/ 

to optimum, detailed comparison between the analytical 
predictions and the experimental measurements have been 
made. The circumferential 'variations of pressure amplitudes 
and phases at several axial positions have been examined in 
straight and variable-area ducts, hard-wall and lined sections 
with and without a mean flow. Reasonable agreement between 
the theoretical and experimental results has been obtained. 
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I. Introduction and Background 

•The prediction of noise emission from aircraft engine- 
duct systems is a problem whose analysis is complicated by 
the nature of the mean flow within the ducts. Such ducts 
carry a high-speed mean flow that possesses streamwise 
gradients due to variations in the duct geometry in addition 
to strong transverse gradients due to the wall boundary 
layers. Although these gradients have been shown to have a 
significant effect on the sound propagation, techniques that 
are capable of predicting their effect in realistic duct 
configurations have not been available until recently, 
primarily due to the difficulty of calculating the effects 
of the streamwise gradients. The reward for the successful 
development of such a technique — that is, the ability to 
identify those duct configurations which possess the most 
favorable sound attenuation characteristics — has prompted 
numerous investigations into the effects of area variations; 
see the review articles^ of Nayfeh, Kaiser, and Telionis 1 and 
Nayfeh?. 

Early studies of wave propagation through variable-area- 
ducts restricted their scope to the case of no mean flow. 

O ' . 

Webster's equation^ for the propagation of a plane wave is 
well-known, and numerous invest igations have verified or 
expanded this theory. Isakovitch^ , Salant 5 and Nayfeh 6,7 

O 

and Nayfeh and Kandil obtained perturbation solutions for 
wave propagation in ducts whose walls vary sinusoidally. 
Several authors (Alfredson^, Lansing and Zorumski "^ 0 and 


1 



2 


11 

Zorumski ) used a discretization technique based on the 

solution for a single duct discontinuity. Nayfeh and 

Telionis applied the method of multiple scales to obtain a 

uniformly-valid perturbation solution for a duct of slowly- 

13 

varying cross section. Beckeme-yer and Eversman used the 
Ritz minimization of fuctionals- with- the governing equations 
as stationary conditions in order to waive the restriction 
of slow variation with axial distance; however, it is not 
known whether their analysis will converge when the duct 
variations are large. Recent studies of the no-mean-flow 
case have been carried out using the method of weighted 
residuals (Eversman, Cook, and Beckemeyer^ ) , integral 
equations (Quinn ) variation of parameters (Kaiser and 

Nayfen 0 ), finite differences (Quinn , Baumeister and 

X 8 IQ 20 

Rice x , and Baumeister y ), and finite elements (Watson and 

21 

Eversman and Astley ) . 

The inclusion of the effects of a mean flow to provide 
a more realistic model makes the problem considerably more 
difficult, and most studies have employed one or more simpli- 
fying assumptions. The most commonly-used assumption is 
that of quasi-one-dimensional flow, which eliminates the 

effect of sound refraction through the boundary layer. 

22 

Powell- studied the propagation of sound discontinuities, 

23 2'4 

while Eisenberg and Kao , Davis and Johnson , Myers and 
25 26 

Callegari , Callegari and Myers , and Nayfeh, .Shaker, and 
27 

Kaiser considered the propagation of the lowest acoustic 
mode through ducts with variable cross sections. Huerre and 
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Karamcheti^ and Grimm^ and King and Karamcheti 3 ^ investi- 
gated the propagation- of the lowest acoustic mode using the 
short wave approximation (ray acoustics) and the method of 
characteristics. Hogge and Ritzi 31 a'ssumed that the duct 
can be broken down into cylindrical and conical sections, 
and they- solved- the wave equation for a uniform mean flow by 
application of the method of finite elements. However, the 
restrictions to one-dimensional flow and/or propagation of 
the fundamental mode limits the usefulness of such analyses 
for realistic situations in which the sound is a combination 
of numerous modes which propagate through a mean flow that 
possesses both axial and transverse velocity gradients. 

Perhaps the only simplifying factor of the problem is 
the fact that, in many practical situations, for example in 
the bypass ducts of high-bypass-ratio jet engines, the duct 
cross-sectional area depends only mildly on the distance 
along the axis of .the duct. Moreover, the growth of the 
mean boundary layer at high Reynolds numbers is a slow 
function of the axial distance, except perhaps at the lead- 
ing edge of the duct inlet. The above factors indicate that 
a perturbation method should be appropriate. SUch a method 
could determine the correction due to weak; nonuniformities 
of the duct shape, a mean flow with a small normal velocity 
component, slow variations of the duct liner properties, and 
a growing boundary-layer thickness, in addition to including 
the effects of large transverse velocity gradients. 



32 

By using the Born approximation, Tam investigated the 
transmission and scattering of spinning acoustic-wave modes 
through a duct nonuniformity. Nayfeh, Telionis, and Lekoudis 
used the method o,f multiple scales in order to study the 
propagation of all acoustic modes in a two-dimensional duct 
with a slowly-varying cross section carrying a sheared 

34 

incompressible mean flow. Nayfeh, Kaiser, and Telionis 
extended the analysis to the more practical annular-duct 
configuration, and Nayfeh and Kaiser 3 ^ applied the technique 
to propagation through a two-dimensional compressible flow. 
These multiple-scales analyses include the effects of stream- 
wise variations of the mean velocity, of the liner proper- 
ties and of the boundary-layer -thickness, and the effect of 
the mean normal velocity and transverse velocity gradient. 
Thus, acoustic propagation through realistic mean flows can 
be examined with the multiple-scale analyses. 

For ducts with larger axial variations, coupling be- 
tween modes will occur. Eversman 36 has developed a proce- 
dure based on the method of weighted residuals to determine 

the reflection effects and modal coupling in : transmission. 

37 

Nayfeh, Shaker, and Kaiser have developed a wave-envelope 
method based on the method of variation of parameters that 
avoids the necessity of having to integrate through each 
axial wave length and is thus more efficient than the stan- 
dard weighted-residual approach. Sigman, Majjigi, and 
Zinn 3 ^ and Abramson 3 ^ developed finite-element models of 
sound propagation in nonuniform ducts carrying compressible 


flows. 
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In contrast to the considerable efforts that have gone 
into the development of analytical and numerical methods to 
predict acoustic propagation in variable-area ducts, few 
experimental studies have been made to test these analytical 
models. Although good agreement between theory and exper— 

• iment has been obtained in lined, straight ducts (for 
example, by Plumblee^ 0 and Plumblee, Dean, Wynne and. Burr in^) , 
no attempt has been made to test the analytical models for 
variable-area ducts in. which the approximations, and assump- 
tions of the models are critical. For this -reason, a pro- 
gram was undertaken in which a primary objective was direct 
comparison of experimental results with, analytical predic- 
tions. The best analytical model available at the start of 
the program was that based -on the method of multiple scales 
and thus this model was chosen for comparison. A second 
objective of the analytical program was -the determination of 
the influence of axial gradients on the "optimum" liner 
properties. As discussed in Section 5.2, under some condi- 
tions, in a variable-area duct the non-parallel results 
exhibit a singularity whose origin has been found to be a 
non— uniformity in the standard parallel— duct theory whenever 
two or more modes have equal eigenvalues. Consequently, the 
perturbation analysis for the variable-area effects are 
not applicable to cases near multiple roots of the parallel- 
duct eigenvalue equation. 
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The experimental facility the instrumentation and the 
design- parameters for the experimental study are described 
in Section 2. In Section 3, an analytical formulation for 
annular, hard-wall ducts with axial variations in the cross- 
sectional area but no mean flow is developed; preliminary 
experimental studies to define conditions for comparison of 
experimental and theoretical results are reported; and the 
detailed comparison of theoretical and experimental results 
is given. The analytical development for the lined test 
section without a mean flow is described in Section 4; in 
addition, the liner properties are determined from the 
experimental data for the constant-area duct, and the 
analytical predictions for the variable-area test section 
are compared with appropriate experimental data. Cases in 
which the duct carries a mean flow are considered in Section 
5; the analytical basis of the computer code for this case 
is reviewed; computational studies of the effects of axial 
variations on wave propagation In ducts with ’'optimum" 
liners are described; and the experimental data in variable- 
area, lined annular ducts are compared with theoretical 
predictions. Section 6 summarizes the investigation. 



2. Experimental Facility and Instrumentation 


2.1 General Approach 

2.1.1 The Measurement Problem 

The basic objective of the experimental investigation 
was to characterize the modal contend of an acoustical 
traveling wave propagating through an annular duct. Changes 
in the acoustic propagation were to be related to changes in 
the geometry of the duct boundaries, the acoustic impedance 
of the duct walls, and the flow of the transmitting fluid 
through the annular passage. 

To allow comparison with theoretical predictions, it 
was necessary to measure both the amplitude and relative 
phase of an acoustical signal at known axial locations at 
the outer wall of the test passage. 

2.1.2 Basic requirements for the experimental system 

The measurement program undertaken necessitated the 
design and construction of a test facility which would allow 
accurate and repeatable acoustic measurements in the pre- 
sence of an appreciable flow. In addition, the test appara- 
tus had to allow for the change of several basic parameters 
with relative ease.* The general design constraints imposed 
on the test-section geometry, air flow, and acoustic-pro- 
pagation parameters will be discussed in turn. 
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The foremost geometric requirement on the test section 
was that it be annular in shape and that it have a varying 
cross-sectional flow area. In addition, it was required 
that at least one of the annular walls be acoustically soft. 

These geometric requirements were met by using a varia- 
ble diameter, acousfcically-hard center-body running through 
an outer cylinder of constant inside diameter. The outer 

wall of one of the annular test sections was made acousti- 

* 

cally soft. An acoustically-hard outer wall and a constant 
area' (uniform) center-body were also constructed. 

Measurements in the presence of flow required the 
establishment of a controlled, constant air flow having a 
flow profile that was reasonably, uniform. Extraneous noise 
was removed from the flow ahead of the test section so that 
an acoustic signal of known frequency, could be detected with 
a good signal- to-noise ratio. -The tests for this_study were 
limited to Mach numbers* less than 0.3 in the variable-area 
test section in order to facilitate comparison with the , 
theoretical results which were 'based on an incompressible 
mean— flow model. However, a capability of handling large 
air flows was desired to allow experimental measurement in 
flow regimes up to a Mach number of 0.6. The air flow . 
requirements were .met by using the compressors and storage 
tanks of a nearby supersonic wind tunnel. The air was 
throttled, muffled, and expanded in a settling chamber 
before passing through a flow straight ener and converging 
section on its way to the test passage. 
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The acoustic requirements dictated that the test sec- 
tion he free from standing waves in the axial direction and 
that provision be made for inserting a pure-tone acoustic 
signal. Furthermore, it was necessary to be able to excite 
various mode shapes at the iHlet of the test sect-ion. The 
minimization of standing waves was accomplished by muffling 
and by the use of horns on both the downstream end of the 
test section and the upstream end of the acoustic signal 
insertion location. . .Mode shape generation was provided by 
using multiple acoustic sources. Each source was controll- 
able in both amplitude and phase, consequently allowing' some 
control of the relative modal amplitudes in the. propagation 
wave . 

The. various aspects of the test apparatus are con-- 
sidere.d in more detail in the following section. ' 



2.2 Experimental Apparatus 
2.2.1 System description. 

The test apparatus is arranged as shown in Fig. 1. Air 
from a high-pressure storage system is regulated to a desir- 
able working pressure by a high-capacity regulating valve 
[1]*'. This commercially available 4-inch pressure-regulating 
valve is designed to reduce an inlet supply pressure of 75 
PSI gage to an exit pressure ranging from 1-30 PSI gage, as 
desired. Under these conditions the valve can deliver 
approximately 10,000 standard CFM of air. These experiments 
were run with a pressure downstream of the regulating valve 
of approximately 25 PSIG. A 4-inch gate valve [2] is used 
to stop and start the air flow. - Final adjustment of the 
fluid velocity is made with a globe valve [3]. 

Two mufflers in series serve to eliminate the flow 
noise and to provide a flow capacitance. Americ an A ir 
Filter ( AAF) Pulsco Model LM18-4 line mute [4] is the first 
muffler in line. This provides, some flow resistance and re- 
•moves the bulk of the throttling noise that results from the 
regulating valves. The flow is. turned and enters an Indus- 
trial Acoustics Company (IAC) model 18-PA-2 silencer [5]. 

This muffler is basically an, 18-inch diameter acoustically 
lined plenum. The diameter was selected to be large enough 
that the flow velocities in this plenum are quite small. It 
thus serves as a settling chamber ahead of the test section. 
*Numbers in brackets refer to items indicated on Fig. 1. 
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Figure' 1. A schematic of test appara 
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The centerbody [6] begins at the upstream end of this muffler. 
It continues at a constant 3-inch diameter through the 
muffler, the contraction section -and the test section. By 
this medns, flow disturbances, which would have resulted 
from the introduction of the centerbody in the established 
flow field near the test section, are avoided. 

After leaving the muffler, the flow passes through a 
straightening device [7] constructed from 8-inch thick 
honeycomb having 1/8 inch 'cells. This device reduces the 
turbulence of the airstream and acts as an additional flow 
resistance to help stabilize the air flow. A contraction 
section [8] then reduces the flow diameter' from about 18 
inches to 6. .100 inches. The contraction section -is basi- 
cally conical in shape, but the downstream end is faired 
into the upstream end of the signal-insertion -section [9] by 
using a long-radius circular arc. 

The signal-insertion section [-9] consists of a 6.100 
inch inside-diameter honed pipe equipped with four loud- 
speaker driver units. These units are equally spaced (90°) 
around the circumference of the outer wall at the same axial 
location. A thin piece of rayon' material covers the driver 
hole opening on the -inner wall of the acoustical-source 
section. The driver units are individually controlled with 
regard to amplitude and phase. The loud-speaker driver, 
units are Atlas Sound Model PD-60, rated at ^2 watts rms with 
a useful frequency range from 70-12K, Hz [10]. The air 
temperature and the flow velocity are measured in the 
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signal-insertion section, which is immediately upstream of 
the test section. 

Two test sections [11] were constructed. These are 
discussed in more detail in Sections 2.2.2 and 2.2.3. 

Two centerbodies [12] were constructed. The straight 
centerbody was constructed from a 3-inch aluminum tubing 
having a 1/4 inch thick wall. A conically-shaped centerbody 
was constructed to provide for the variable-flow area. The 
details of the construction of this centerbody are given in 
Section 2.2.4. 

A reference microphone probe [13] located at the en- 
trance of the test section is used to maintain the sound 
pressure level (SPL) constant at the reference point. The 
signal from this probe also serves as a reference for deter- 
mining the phase of the signal from the measuring probes. 

The measuring probes [14], which are flush mounted, pro- 
vide the acoustic pressure at the outer wall of the test 
section. Six axial probe positions are located in each test 
section. Sufficient circumferential positions are available 
at the entrance of the test section to define the modal 
content up to the third circumferential mode. The details 
of the probe locations are explained in Section 2. 2.. 2, which 
defines the geometry and construction of the test sections. 

The test section is terminated with a horn [15] which 
was fabricated from the outer section of a large outdoor 
speaker and, therefore, is not ideally exponential. The 
horn does have a high efficiency at the desired test 
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frequencies so that no standing waves were detected in the 

test section (see Section 3-2). 

The system exhausts into an anechoic chamber [l6]. The 
anechoic chamber is of plywood construction and is lined 
with a 4-inch reticulated foam. Provision is made to ex- 
haust the' air through the top of the chamber. All lines of 
sight into the nozzle exit are blocked. The system is 
anechoic to 125 Hz. Experiments indicated that the addition 
of the anechoic termination greatly reduced upstream sound 
propagation into the test section from external sources. ■ 

2.2.2 Hard-walled test section 

An acoustically-hard test section was fabricated from a 
standard 6-inch schedule 40 pipe. The pipe was bored to. the 
required 6.100-inch inside diameter. In addition,, the test 
section was carefully honed to insure that the required 
diameter is constant around the circumference and along the 

axial length. 

The. geometry of the variable-area centerbody ' dictates 
the axial positions for the microphone probe placements. 
Determination of the modal content of the propagating wave . 
dictates the pircumf erential- probe placements. Six axial 
probe stations are provided, as shown in Fig. 2. The number 
of circumferential probe locations vary, depending on the 
axial location. 

Station 1 (at the upstream end of the test section) 
serves as the reference station. It has eight circumferen— 
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tlal probe positons. When a variable-area centerbody Is 
used, it is positioned so that station 1 is just upstream of 
the increasing-diameter portion of the centerbody. Station 
2 has two diametrically-opposed circumferential probe posi- 
tions and is axially located approximately midway between 
the minimum and maximum diameter locations on the variable- 
area centerbody. Station 3 has four equally-spaced circum- 
ferential probe positions. Its 1 axial location corresponds 
to a position just prior to the location of the maximum 
diameter of the variable-area centerbody. 

All probe holes in the hard-wall section were machined 
such that one probe assembly could be moved from location 
to location and the probe tip would be flush with the inner 
wall to within an acceptable tolerance. 

2.2.3 Soft-walled test section 

The computer code described in Section 5*1 was used to 

examine the effects of several liner constructions on the 

42 

acoustic propagation as a preliminary step in the design . - 
The influence of the liner was investigated for both straight 
and. variable-area centerbodies with and without a mean flow. 
The liner design was tuned to approximately 3000 Hz, but the 
peak attenuation was not made so large as to reduce the 
input signal below the flow noise in the test section. The 
selected soft-walled test section consists of a hard outer 
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shell fabricated from an 8-inch diameter steel tube and a 
0 . 625-inch thick "Flex-Core" honeycomb which separates the 
facing sheet from the outer shell (see Fig. 3). The porous 
facing sheet is made from "Feltmetal", a material supplied 
by the Brunswick Corporation. The material selected has a 
nominal resistance of 40 CGS Rayls. However, visual inspec- 
tion indicated some obvious variations in porosity across 
the surface of this sheet. This variation of the facing- 
sheet resistance is a major difficulty in obtaining agree- 
ment between the results of the experimental testing program 
and the analytical predictions for the soft-walled test 
cases . 

One of the major construction difficulties faced in 
this project was that of bonding the Feltmetal to the honey- 
comb without destroying the porosity of the Feltmetal. The 
construction technique finally worked out to preserve the 
porosity of the sheet is as follows. A sheet of thermal 
adhesive is placed on the flex-core. Then a small hole is 
punched in the sheet adhesive over each flex-core cell. The 
adhesive is next heated by a heat lamp, which causes the 
adhesive to contract on the lip of each cell. The test 
section is then assembled and cured under pressure. Exami- 
nation of the final assembly has shown that at least 95 % of 
the lips of individual cells are bonded to the Feltmetal by 
this technique without appreciably destroying the porosity 
of the Feltmetal. 
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The axial placement of probes is dictated by the 
variable-area centerbody and is identical to that employed 
in the hard-walled test section as discussed above. How- 

i 

ever a early tests with the hard-walled test section indi- 
cated that more circumferential probe positions at stations 
one, three and six would be desirable. Further, these probe 
positions should not be equally spaced around the circum- 
ference. Consequently axial stations one and three of the 
soft-walled test section have ten (10) probe positions as 
specified in Fig. 3. 

In order to assure proper alignment of the. probe tips 
flush with the Feltmetal liner surface in the soft-walled 
experiments, separate probe tips are permanently left in 
each probe location, and the microphone and probe block are 
moved from tip to tip during the course of the experiment. 
Great care has been taken to insure that the probe tips are 
identical in size and have similar acoustical characteris- 
tics. Each probe tip was calibrated with respect to both 
amplitude and relative phase prior to installation into the 
soft-walled test section. Section 2.3.2 and Figure 6 pro- 
vide more detail about the probe assembly. 

2.2.4 Variable-area centerbody 

The variable-area centerbody was designed to produce a 
significant velocity gradient along its axial length. The 
computer code described in Section 5*1 was used to examine 
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the effects of the axial gradients on the acoustic propaga- 
tion through the test .section. The area variation was 
selected large enough to produce a measurable non-parallel 
effect on the acoustic disturbance but not so large as to 
invalidate the theoretical results with which experimental 
data were to be compared. The essential dimensions of the 
variable— area centerbody that was selected are shown in Pig. 
4. The approach diameter is two' inches, and the diameter 
increases to a 4.296-inch maximum. The variable-area cen- 
terbody tapers to. the 3-inch exit diameter at approximately 
half the approach rate. The various mating surfaces have 
been blended, and the entire' surface is polished to minimize 
acoustic or flow perturbations that might result from' sur- 
face imperfections. 

For the soft-walled tests, the variable-area centerbody 
was positioned so that the start of the incline was just' 
downstream from the test-section probe reference station 
(station 1). The maximum centerbody diameter was then' 
slightly downstream from the third axial probe station 
(station 3) s so that measurements were made prior to the 
throat . 

The exact centerbody geometry was carefully measured 
with a dial indicator, and. the- results are shown in Fig. 4. 
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2.3 Instrumentation 


The instrumentation associated with the experimental 
facility consists of three parts:' source instrumentation, 
servo-control instrumentation and probe instrumentation. 

Each of the three parts is explained in detail. A schematic 
of the complete electronic system is shown in Fig. 5- Since 
much of the instrumentation has evolved from the Spectral 
•Dynamics vibrational control system, this system is explain- 
ed in detail prior to explaining each of the separate in- 
strumentation sections. 

2.3.1 Basic control system 

The basic instrumentation system has been built around 
a system of Spectral Dynamics instruments that are designed 
to control vibrational shakers. The four instruments are an 
SD 104A-5 sweep oscillator, an SD 105-A amplitude servo /monitor , 
an SD 122 dual channel tracking filter, and an SD 127 MZ/TFA 
control. A brief explanation of the function of each of the 
instruments and their interconnection follows. 

The SD 104 a- 5 sweep oscillator provides a constant 
amplitude sinusoidal voltage. This voltage is used as a’ 
source for the SD 105-A amplitude servo/moriitor and as a 
reference frequency by the SD 122 tracking filter. 

The SD .105-A amplitude' servo/monitor is a servo-control 
device designed to maintain either acceleration, velocity or 
displacement constant in shaker 'application. No instrument 
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modification is required for use in these acoustical experi- 
ments. In these experiments the SD 105-A is used to keep 
the' sound pressure level constant at the 12:00 o'clock 
position of axial station 1. A "cross-over" feature is also 
utilized to insure that the power rating of the acoustical 
drivers is not exceeded. 

The SD 122 tracking filter is a dual channel tunable 
filter. The unit uses the constant amplitude output of the 
SD 104A-5 sweep oscillator as a reference frequency. Band- 
widths of 2, 5, 10, 20, 50, and 100 Hz are available. In . 
these experiments, one channel of the tracking filter is used 
in the serro control circuit (10 Hz bandwidth), and the other 
channel is used in the probe measurement circuit (2 Hz band- 
width) to separate the imput frequency from the flow noise. 

The SD 127 MZ/TF-A control is basically an integrator 
used to obtain impedance plots in vibration applications. 

The SD 127 is interconnected with the SD 122 tracking filter 
and serves as a phase meter in the experiments .- 

2.3.2 Microphone Probe .Assembly 

One of the crucial parts of the measurement system is 
the microphone probe assembly. A schematic of the assembly 
is shown in Fig. 6. The four basic components of the prpbe 
assembly are a 1/2-inch B&K type 4133 microphone, a probe 
block, a termination coil,, and a probe tip. The basic 




Figure 6. Microphone probe assembly 
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diameter of the probe tip is 1/8-inch. The design is based 

on a similar probe utilized by General Electric in jet 

43 

engine studies 

Although the purpose of the termination coil in this 

design was to eliminate resonances in the microphone assembly, 

resonances similar to those observed with the standard B & K 
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microphone probe assembly were present. The response of 
such a probe may be improved over a narrow frequency band by 
suitable damping. However, since only three discrete fre- 
quencies have been examined in the experiments., the probe 
simply ha,s been calibrated at the frequencies of interest, 
rather than attempting to use damping. 

The microphone probe assemblies have been calibrated 
using a standard B&K microphone probe calibrator, using 
standard B&K procedures as follows. Care is taken to. insure 
that .the pro.be tip is maintained in the uniform sound field 
produced by the calibrator. With the probe assembly in- 
serted into the calibration assembly, the frequency of the 
sine random generator is set using a counter. The output of 
the B&K type 1024 sine random generator Is fully compressed 
and compression slowly reduced until the desired calibration 
level is reached. The rms voltage of the reference micro- 
phone and the probe assembly are recorded. The effective 
sensitivity of each of the microphone probe assemblies is 
calculated from the following formula: 
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where S p = Sensitivity of probe microphone (N/M 2 ) 

= Sensitivity of reference microphone (N/M 2 ) 

V r = RMS voltage of reference microphone 
Vp = RMS voltage of probe microphone 

Note that since the same amplifier is used in measuring both 
the reference microphone and probe microphone voltages, any 
errors associated with the amplifier gain will cancel. In 
these experiments, the microphone probe assemblies have been 
calibrated at 840, 26*10 and 3150 Hz. 

In the hard-walled experiments it was possible to use 
one microphone probe assembly. This reduced the amplitude 
and phase errors associated with using many probes and 
amplifiers. In the experiments using a soft-walled test 
section one microphone probe block has been used together 
with many permanently-installed microphone probe tips. 

Care has been taken to insure that the length of each probe 
tip is within a + .010-inch tolerance. Each probe tip has 
been individually calibrated in the probe block with respect 
to both amplitude and phase. The sensitivity of each of the 
probe tips is within 5 % of the average value, and the phase 
shift associated with individual probe tips is' within 
the + 2 degrees recording accuracy. 



2.3-3 Acoustical source electronics 


.As described previously , the electronic driver instru- 
mentation consists of the SD lO^A-5 sweep oscillator, the 
SD 104-A amplitude servor-monitor , one channel of the SD 122 
tracking filter, a phase shifter, a power amplifier, the 
reference microphone probe and a signal conditioning amplifier. 

The output of the reference probe-microphone assembly 
is fed via the SD 122 tracking filter (with a 10 Hz band- 
width) to the servo-input of the SD '105-A servo-monitor. An 
initial voltage setting is -then selected which is within 
the rating of the acoustical 'drivers. The servo-control 
then maintains the SPL constant at the reference microphone 
probe . 

2.3-^ Data acquisition system 

The data acquisition- electronics consist of the data 
microphone probe assembly and its signal conditioning, ampli- 
fier, the SD 122 tracking filter and the SD 127 controller. 

The microphone probe signal is filtered and amplified by the 
signal conditioning amplifier. Limited filtering is effected 
by this amplifier. The filtering for data acquisition is 
accomplished by use of the SD 122 tracking filter. A band- 
width of 2 Hz has been used in the SD 122 tracking filter. 

Thus, data has been obtained only within + 1 Hz of the 
signal supplied to the acoustic drivers. 



The rms value of the amplified and conditioned output 
is measured with a true rms voltmeter. Under no-flow 
conditions relative rms values of amplitude are limited only 
by the- accuracy of the' digital voltmeter. This error is 
typically less than . 1 %. 

Phase, relative to the reference probe, is obtained by 
utilizing the phase capabilities of the SD 127 MZ/TFA con- 
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troller and can be recorded to an accuracy of + 2° . 



2.4 Experimental Set-up 


In some of the experimental tests an attempt was made 
to minimize the presence of higher modes in the acoustic 
signal that enters the test section. The oscillator was set 
to the chosen frequency, which was verified with the aid of 
a counter, and the servo controls that regulate the sound 
pressure level amplitude, were by-passed. Each of the four 
acoustic driver units was tested independently to determine 
the amplifier and phase-shifter settings required to produce 
specified signals at a measuring station. The measuring 
station for this set-up procedure was at the same circum- 
ferential location as the acoustic driver being investigated 
For the soft-walled test' section, these set-up measurements 
were made at the downstream- end of the signal-insert ion 
section, and for the the hard-walled test section the 
measurements were made -at the first station- (station 1) of 
the test section. After the capabilities of each acoustic 
driver unit had been determined, each was supplied, in turn, 
with the same signal' from the oscillator, and the amplifier 
and phase shifter for each driver unit were set to produce 
identical acoustic signals at each unit’s measuring station. 

All four acoustic drivers were then connected in 
parallel, to be driven by the same oscillator, and the 
system returned to servo control following the normal pro- 
cedures outlined in the SD manual. The servo control 
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operated to maintain the sound pressure lev.el constant 
during the tests. 

. There was, of course, some interaction between the 
acoustic drivers. The impedance seen by one acoustic driver 
was altered by the' ac.oustic signals coming from the other 
three acoustic drivers. -Consequently, the phase of the 
signal obtained was somewhat different from -the arbitrarily 
selected phase angle used in the set-up procedure. Thus, 
the dominant acoustic signal was the plane mode, with 
relatively smaller amplitudes of the spinning modes present. 



2 . 5 Summary 


An acoustic duct facility has been constructed for 
studies of sound propagation through mean flows that possess 
both axial and transverse velocity gradients. The design of 
the duct facility controls the mean flow so that the mean 
boundary layer in the test section is thin, and controls the 
ambient noise levels through the use of line mufflers and an 
anechoic termination. The annular test section can be 
selected to have either a rigid or lined outer wall; the 
cross-sectional area can be either constant or variable by 
selection of a straight or variable-radius inner wall, which 
is rigid. The duct liner consists of a 40-Rayl Peltmetal 
facing sheet that covers .625-inch cellular cavities. The 
duct has been instrumented to permit acquisition of narrow 
band data for determination of the relative amplitudes and 
phases of the acoustic disturbance along the outer wall of 
the test section. The test sections are provided with 
sufficient microphone probe locations to permit resolution 
of the modal content of the acoustic signal up to the third 
circumferential mode. 
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3. Hard-Walled Annular Ducts Without Plow 


In this section, an analytic solution is determined to 
the problem of sound propagation in hard-walled annular 
ducts with slowly-varying cross sectors in the absence of 
mean flow. Experimental studies of propagation in the hard- 
walled duct are described, and the experimental and analyti- 
cal- results are compared. 

3.1 Theory 


In the absence of mean flow and for negligible dissi- 
pation, the propagation of sound in a hard-walled annular 
duct can be mathematically described by the wave equation 

( 1 ) 

subject to the boundary conditions 


1 3^<f> 

c^ 3t ' 


V 2 <f>- = 0 


!r = eRI fx at r = Rl ^ Xl ) 


3<j> 

3x 

3<j) 


= eR£ g at r = R 2 (xi) 


(2) 

' (3) 


where <j>(x,r, 0 ,t) is the potential function, e is a small 

dimensionless parameter, xi = ex is a long scale describing 

the slow axial variations of the duct cross-section, Ri and R 2 

are the inner and outer radii of the duct, and primes denote 

differentiation with respect to the argument. 

An approximate solution to Eqs. (l)-(3) is sought by 

ii 5 

using the method of multiple scales - in the form 

<J>(x,r,e,t) = [ipo(xi,r) + eih(x 1} r) + ...]exp[i(ri +m0)] 

( 4 ) 
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(5) 


3{t 

The spinning mode number m is an integer and 
W'~ ~ w s 3x = k ^ Xl ^ 

where to is the sound frequency, which is. assumed to be 
constant, and k(xi) is the wavenumber. The axial and time 
derivatives become 

■3_ = a_ 

3t 3n (6) 

a_ „ k 3_ + 3_ 

3x 3r) 3xi 

Substituting Eqs.. (4)- (6) into Eqs. (l)-(3) and equating 
coefficients of like powers of e, one obtains 
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= A mn (xi) Vn (r ’ Xl) 

J 1 ( T Ri) 

^mn ~ J m < ' Y mn r ' ) “ Y m^ mn R i ) Y m^ Y mn r ^ 


The y are the roots of 
mn 


J m (YR2)Y m (YRl) " J m (YRl)Y m (YR2 ^ = 0 


(13) 

(14) 


(15) 


and the wavenumber k is- related to y bv 

mn 'mn J 



( 16 ) 


In Eq. (13), A mn ( x i) is undetermined at this level of appro- 
ximation; it is determined by imposing a solvability condi- 
tion at the next level of approximation. 

With ip 0 and. k mn known, the inhomogeneous Eqs. (10)-(12) 
have a -solution if , and only if, their right-hand sides 
satisfy a solvability condition, since their homogeneous 
parts have a nontrivial solution. To determine this solva- 
bility condition, multiply Eq. (10) by rifi 0 ,, integrate the 
result by parts from r = R x to r = R 2 to transfer the deri- 


vatives 'from to o , and use Eqs. ( 7 ) - ( 9 ) to obtain 



- (k mn i^) dr . - . .(17) 
Substituting for 9<h/3r at r = Ri and R 2 from Eqs. (11) and 
(12) into Eq. (17), we have 
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Using Leibnitz* rule, we rewrite Eq. (18) as 
R2 


1x7 [k mn = 0 


(19) 


or 
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Since 9i|/o/3r = 0 at r = Ri and R 2 , 
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Hence, 
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where xi r is a reference axial location, the complex ampli- 
tudes A (x, ) at the reference location must be found 
experimentally, and 
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Thus, the complex acoustic pressure variation of a single 
mode at the outer wall is given by 
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P mn U,R 2 ,t) = 
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where the eigenfunction: ¥ is defined by Eq. (1*0 arid the 
axial amplitude variation is defined by Eq. (22). This may 
be re-written as 


„ ^ „ D nrn (lt 'r> 
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where the constant is determined experimentally and is 
the complex amplitude of the. mode at the reference station. 
The acoustic pressure is summed over all modes present in 
the duct, and the result is expressed in terms of an ampli- 
tude and phase of the total signal for comparison with the 
experimental data in Section 3.3. 


exp[i( 


/ 



3.2 Preliminary Investigations - Using the Hard-Walled Test 

Section 

To facilitate the comparison of experimental and 
theoretical results, it is desirable to obtain experimental 
data for conditions which .are not close to a mode cut-on ' 
frequency. By considering only well cut-on modes, strong 
, reflections leading to upstream propagating modes and strong 
axial standing waves are avoided. Thus a preliminary investi 
gation was conducted to determine the conditions for which 
strong standing waves occur for cases with and without mean 
flow. These results are easily compared with the theoretical 
hard-walled cut-on frequencies, giving a relatively simple 
check of the instrumentation and the validity of the theore- 
tical modelling of the experimental facility. Finally, the 
frequencies at which more extensive and detailed tests were 
to be conducted are identified in this section. 

3-2.1 Low- frequency standing waves 

Since the acoustical horn on the end of the duct is 
least efficient at low frequencies, the system has been 
checked for .standing waves at a frequency below the cut- 
on frequency of the first circumferential mode. .A frequency 
of 840 Hz was selected for 'this examination, and the ex- 
periments were conducted using the straight centerbody. 
Acoustic pressures were measured at each axial location and 
were found to be within + 0.5 dB in amplitude. It was 
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concluded that no axial standing waves existed. Since the 
efficiency of the acoustical horn would increase with in- 
creasing frequency, it was also concluded, that no substan- 
tial standing waves existed in the frequency range of in- 
terest, provided that all propagating modes were well cut 
on. 


3.2.2 Cut On of Higher .Circumferential Modes 

At frequencies close to the cut-on frequencies of the 
circumferential modes, strong reflections from the duct 
termination are anticipated, and it was desired to avoid 
conducting the experiments at these frequencies. 'The cut-on 
frequency is- defined as the frequency at which the axial 
wavenumber for that mode is zero; that is, the wave is not a 
propagating wave. For a circular, annular duct having hard 
walls, the acoustic waves at cut on may be visualized .as m- 
lpbed circumferential spinning waves of constant amplitude 
along the axial length. 

ill 

Plumblee et al have shown that the cut-on frequencies 
may be obtained experimentally by sweeping an acoustic . 
driver through the appropriate frequency range and monitor- 
ing the pressure amplitude at some downstream station. 

Using this procedure, the strong reflections that occur near 
cut-on frequencies generate large pressure amplitudes within 
the duct which are easily detected. 



Figure 7 gives the voltage output from the measuring- 
microphone amplifier, and thus provides an indication of the 
pressure amplitude as a function of the exciting frequency. 
These measurements were made in the hard-walled test section 
with the straight centerbody (1.5-inch radius) installed. 

The cut-on frequencies correspond to the peaks and are 
clearly visible. The amplitudes shown in Fig. 7 have not 
been corrected for the response of the microphone probe. 
Since the relation between the voltage output and the pres- 
sure amplitude is not constant over the entire frequency 
range, the results in Fig. 7 cannot be' used in a quantita- 
tive manner to calculate sound pressure amplitudes. All . 
the peaks observed in Fig. 7 agree with those expected for 
'the hard-walled annular duct. Theoretically, the cut-on 
frequencies may be obtained by setting the wavenumber, k^, 
equal to zero in the eigenvalue equation 


k 2 = SL. - y 2 
mn c ■ 'mn 


(25) 


where to is the circular frequency, c is the speed of sound 
of the medium, and y mn - is the characteristic number. The 
characteristic numbers, Y mn » are the roots of the .equation 


t 

J^CyROy^yRz) - J i ;(yR 2 )y i ;(yRi) = o 


( 26 ) 


The freqencies satisfying' this equation (the cut-on fre- 
quencies) for the hard-walled test section are given in 


Table 1. 



Voltage, output 



Figure 7. Frequency characterization of hard-wall test section with uniform centerbodv 

installed. 
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Table 1. Theoretical' cut-on frequencies for a hard-walled 
duct with R 2 = 3.05". 


Mode 

Number* 


m,n 


Cut-on Frequencies, Hz 


Inner radius 
R i" = 1.0” . 


Inner radius 
Ri = 1.5" 


0 , 0 
0,1 
0 , 2 - 
0,3 
0,4 


3356 

6518 

9713 

12920 


- 4336 
8558 
12802 


1,0 

1,1 

1,2 

1.3 

1.4 


1066 

3615 

6639. 

9791 

12977 


2,0 

2,1 

2,2 

2,3 


2024 

4317 

7002 

10024 


3.0 

3.1 

3.2 

3.3 


2870 

5262 

75.97 

10412 


939 
4459 
' 8616 
12840 


1856 

4813 

8789 

12955 

2736 

5366 

9073 

13143 


4.0 

4.1 

4.2 

4.3 


3657 

6252 

8405 

10954 


3573 

6071 

9461 

13405 


5.0 

5.1 

5.2 

5.3 


4417 

7185 

9356 

11651 


4372 

6880 

9946 

13755 


6 , 0 
6,1 
6,2 


5167 

8061 

1,0345 


5143 

7748 

10521 


*m - Circumferential mode number 
n - Radial mode number 
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The agreement between the theoretical and experimental 
values is good. The resonant frequencies, which are- appro- 
ximately three percent higher than the theoretical cut-on 
frequencies, are sharply def-ined. The three percent shift 
in the frequency of peak response is a consequence of damp- 
ing and is . consistent with the observations of an earlier 
study (Ref. 4l, pp. 53-55). For this particular hub-to-tip 
ratio, no- radial modes were encountered before the third 
circumferential mode cut-on frequency. 

Figure 8 is a similar plot showing the variation of the 
amplitude with frequency for the system with the variable- 
area centerbody installed. In the signal-insertion section, 
the variable centerbody has a radius of 1 inch. This smaller 
hub-to-tip ratio shifts the cut-on frequencies of the funda- 
mental radial modes (n=0) to higher values and lowers the 
cut-on frequencies of the higher radial modes (see Table 1). 
The shift of the resonance of the spinning modes to higher 
frequencies is apparent in a comparison o,f Figs. 7 and 8. 

As before, the resonance peaks appear at frequencies that 
are three percent higher than the hard-walled cut-on fre- 
quencies. The first radial-mode cut-on frequency in the 
signal-insertion section is lowered from. 4336 Hz to 3356 Hz, 
and this appears as a sharp valley in the amplitude spectrum 
at 3450 Hz (three percent higher than 3356 Hz). Although 
the first radial mode is cut on in the. signal— insertion sec- 
tion at this frequency, it cuts off in the variable-area 
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Frequency Hz 

Fiqure 8. Frequency characterization of hard-wall -test section with variable 
area centerbody installed. 
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section just downstream; the reflection from the variable- 
area centerbody at the cut-off point apparently causes .a 
partial cancellation of the signal upstream. The sharp 
valley at 3720 Hz corresponds to cut on of the (1,1) mode. 
The variable-area centerbody introduces a slight resonance 
in the region of 2300^2450 Hz; the origin of this resonance 
is unclear, and it was decided to avoid further testing in 
this range. Using the results of these two figures and 
attempting to avoid resonance conditions, we decided to 
carry out detailed experimental studies at three frequen- 
cies: ' one, 840 Hz, just below cut on of .the first spinning 
mode; a second, 2640 Hz, just before cut on of the third 
spinning mode; a third, 3150 Hz, between cut on -of the third 
spinning mode and the first radial mode. 

3.2.3. Preliminary Tests with Mean Flow 

When air is flowing through the test section, the cut- 
off frequencies are changed. Theory predicts that the 
effective cut-on frequency will vary with flow velocity as 

f ’ = V 1 - M 2 (27) 

where f is the original eut-on frequency without flow. 

Tests were run using the variable-area centerbody in 
both the soft-walled and hard-walled test sections. The 
tests confirmed the shifts in the resonance frequencies and 
showed a strong resonance for the first radial mode, unlike 



the no-mean-flow cases. The test frequency of 3150 Hz was 
selected in an attempt to avoid this strong resonance. 

Measurements also were made to determine the effect of 
the mean flow on the control of the input signals. As noted 
in Section 2, the input drivers could be controlled to 
minimize the presence of .the circumferential modes; the 
acoustical wave at station 1 could be maintained, plane 
within .+ 2 dB in amplitude and + 15° in phase for the 3150 
Hz test frequency, and- within + 0.1 dB in amplitude and a 
+ 2° in phase f.or the 840 Hz signal. Because of the dif- 
ference in impedance presented to the sound sources, the 
relative amplitudes of the higher modes increased somewhat 
as the free-stream velocity increased. 



3.3. Comparison of Experimental and Th'eoretical Results 

The theoretical prediction method described in Section 
3.1 predicts the attenuation and phase shift of each acous- 
tic mode as it propagates through the test section, provided 
the , relative amplitude of each mode is known. 

In order to establish the modal content of the signal, 
detailed measurements of amplitude and phase around the 
circumference at station 1 are necessary. Since the highest 
test frequency (3150 Hz) allows the third circumferential 

I 

mode to propagate, at least seven measurements of amplitude 
•and phase are required. The duct instrumentation .was de- 
signed -to provide this data. A least.-squares technique- is 
used to obtain the complex constants B mn of Eq. (24). Once 
these constants are determined, the acoustical pressure at 
the outer wall at the. downstream locations are obtained 
using Eq. (24 ) . 

.The first experiments performed .were those using a 

hard-walled test section and a variable-area centerbody. 

Figures 9a and 9b show the circumferential variation of the 

pressure - amplitude and phase, respectively, at the outer 

walls for three axial stations for an input frequency of 

840 Hz. 'Only one mode, the plane mode, is cut on at this 

frequency. Both the amplitude and the phase measurements 

* 

are in good agreement with the theoretical predictions. The 
theory shows the intensification of the signal as it pro- 
pagates from station 1 to station 3 although it' slightly 
overpredicts the increase. 



Pressure Amplitude, N/M 



Figure 9a. Comparison of the theoretical and experimental pressure amplitudes for 
the hard-walled case; M = 0 and f - 840Hz 
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Figure 9b. Comparison of the theoretical and experimental pressure phases for the 
hard-walled case; M = 0 and f = 840 Hz. 
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Figures 10a and 10t> are similar plots for an input 
frequency of 2640 Hz. Three modes (plane, 1st circumfer- 
ential, and 2nd circumferential) are cut on at this fre- 
quency. Both the amplitude and the phase predictions agree 
very well with the measured values. The measured increase 
in the acoustic intensity as well as the decrease in the 
circumferential variation of the pressure are both predicted 
by the theory. 

Figures 11a and 11b show results for an input frequency 
of 3150 Hz. Four modes (plane, 1st, 2nd, and 3rd circum- 
ferential) are cut on at this frequency. Although there is 
a reasonable agreement between the measured and predicted 
phases., the circumferential variation of the pressure ampli- 
tudes is not adequately predicted by the theory. Neverthe- 
less, the theory does correctly show that the circumfer- 
ential average of the pressure amplitude increases as the 
signal propagates down the duct. 

It is apparent from all three cases that more probe 
positions are needed at downstream locations especially if 
significant components of higher spinning modes are present. 
Consequently, the lined-duct tests, which were conducted 
later, included more probe locations at stations 2 and 3 . 
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Figure 10a. Comparison of the theoretical and experimental pressure amplitudes for 
the hard-walled case; M = 0 and f = 2640 Hz. 
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Comparison of the theoretical and experimental pressure phases for 
the hard-walled case; M = 0 and f = 2640 Hz. 
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•Figure 11a. Comparison of the theoretical and experimental pressure amplitudes 
- for the hard-walled case; M = 0 and f = 3150 Hz. 
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Figure lib. Comparison of the theoretical and experimental pressure phases for the 
hard-walled case; M = 0 and f = 3150 Hz. 


4.. Lined Annular Ducts Without Flow 


In this section, an analytic expression is derived for 
the propagation and attenuation of sound waves in annular 
lined ducts having slowly-varying cross sections and in the 
absence of mean flow. Experimental studies of propagation 
in the lined duct are described, and the experimental and 
theoretical results are compared. 

4..1 Theoretical Considerations 

In the absence of a mean flow and viscous dissipation, 
the mathematical statement' of the problem can- be written as 


n 2 j. 1 

V + - S* 3# - 0 




(28) 

3 4> Bi (xi ) 34) _ 

3r c 3t 

eRJ 

9x 

at r =’ Ri (xi ) 

( 29 ) 

3<j) . 62 (Xi ) dip _ 

3r c 3t 

eRi 

3x 

at r = R 2 (x 1 ) 

( 30 ) 


where ^(xj) and 3 2 (xi) are the admittances of the two 
walls. 

An approximate solution is sought to Eqs. (28 ) — (30) by 
using the method of multiple scales in the form 

4>(x,r,0,t) = [iJj o (xi ,r ) + eipi(xi,r) + ...]exp[i(n +m0)] 

(3D 

where m is an integer and 

|a = - |n = k(Xl) (32) 
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with w being the sound frequency; k is the complex wave- 
number in the streamwise direction whose imaginary part 
represents the attenuation of the sound. Substituting Eqs. 
(31) and (32) into Eqs. (28)-C30), using Eq. (6), and equating 
coefficients of like powers- of e, we obtain 


Order e 0 

9 2 i|j o , 1 

r 3r 

|fe + 



= 0 at r = Ri 
= 0 at r = R 2 


0 


( 33 ) 

( 34 ) 


( 35 -) 


Order e 


8 2 <Ki . , 1 9^i 
+ r 3r 


+ 



k 2 ; 



2ik 


o 

9xi 


ik'i^o 

( 36 ) 


§f ~ + = ikRIi^o- at r = Ri . (37)' 

||j. _ igJLz. ^ = ikR^ 0 at r = R 2 ( 38 ) 


The solution of Eqs. (3 3 ) — ( 3 5 ) is 

1^0 = A(x! ) Y(r ;xi ) (39) 


where 


y = J (v r) 
m f mn ' 


Y J* (y Ri ) + 
_mn m ’ran 1 


i(d$ 


L J (Y_Ri) 


Y Y 1 (y Ri ) + 
'mn m 'mn 1 ' 


c m mn Y (v r) 

Bn 


nr 'mn 



(40) 

( 41 ) 
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and the y are the roots o,f 
mn 

[yJ^(yR 2 ) - J m (YR 2 )][yY^(yR 1 ) + YjyRi)] 

- C.yJ^CyRi) + J m (YRi')][YY;(YR?'i 

-^• Y m ( YR 2 ) ] . (42) 


With bo and k mn known, Eqs. ( 36 ) — ( 38 ) have a solution 
if, and only if, their inhomogeneous parts satisfy a solva- 
bility condition: To determine this solvability condition, 

multiply Eq. ( 36 ) by rbo, integrate the result by parts from 
r “ Ri to r = R 2 to transfer the derivatives from to bo, 
use Eq. (33), and obtain 


Obo 


iii 

9r 


rbi 


3 b 0 1 
9r J 


r=R 2 

r=Ri 


R 2 

i r (kbo)dr 

Ri 


(43) 


Substitution for and,-||—- at r = Ri and R 2 from Eqs. 

(34), (35), (37), and ( 38 ) into Eq.- (43) results in 


R; 


k R 2 RJbo 

mn Y 


r=Rj 


k RiRlbo 

mn 1 1 Y u 


r=Ri 


Ri 


9xi 


( k b 0 ) dr = 
mn Y 

(44) 


Using Leibnitz’ rule and the fact that bo 
Eq. (44) as 


9 

9xi 


Hence, 


R 2 

[A mn k mn rT2dr] 

Ri 


0 


A rewrite 

m n * 


(45) 
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A 


mn 


D (xi ) 

p (x — fl nrn (x ’r> 
mn 


(46) 


where xi r is a reference station, the complex amplitudes 
A (xi ) at the reference location are determined experi- 
mentally, and 


D. 


2 _ 
'mn 


k mn R| 


2y‘ 


Ct 


co 2 3 1 


mn 


mn 


m 

c^ " Rl 


k R? 
mn 


]T 2 ( r 2 ;xj ) - ^ 

■ m 


mn 


[y 2 

*■ mn 


w 2 3i 


Ri 


■]¥ 2 (Ri.;x 1 ) 


(47) 


Analogous to the discussion in Section 3> the development of 
the acoustic pressure at the outer wall leads to 


P ™ lX,r2 ’ ’ “ V (x ‘> V (R ^ Xl r > 

X 


exp[i ( 


k dx + m0 - tot)] 
mn 


(48) 


where B is the complex mode amplitude at the reference 
mn 

station. The pressure is summed over all modes, and the 
result is expressed in terms of an amplitude and phase for 
the total signal. These amplitudes and phases are compared 
with the experimental data in Section 4.3. 



4.2 Experimental Investigation 
4.2.1 Introduction 

The soft-walled experiments were performed using the 
apparatus as described in Section 2. Both straight and 
tapered centerbodies were used so that cases for both vari- 
able area and constant area were included. 

The introduction of a soft-walled test section brings 
with it several complications. First, the impedance of the 
wall. is no longer well-known. Also, the insertion of a 
soft-walled section in a tube which is otherwise hard-walled 
gives rise to impedance discontinuites (and thus reflec- 
tions) at the hard-wall to soft-wall joints. 

Thus, the three necessary steps in this investigation 
are (1) determine if any problems with standing waves 
occur, (2) find the impedance of the soft-walled section as 
installed and calculate the basic facing-sheet properties 
(resistance and characteristic frequency) from the imped- 
ance, and (3) measure the -acoustic pressures and phases for 
the. variable-area case. The results from this last step are 
compared with the theoretical predictions. 

4.2.2 Standing-Wave Investigation 

The experimental approach to standing-wave detection is 
basically the same as that outlined for the hard-walled 
experiments. An acoustical driver, located at the 12:00 
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o'clock position, is operated at constant voltage as the 
frequency is slowly swept from 500 Hz to 5000 Hz. An 
acoustical probe simultaneously detects the sound-pressure 
level at the discrete axial locations. The amplified micro- 
phone output is passed through a 10 Hz tracking filter and 
it is recorded. The result is an uncalibrated pressure 
amplitude vs. frequency plot with the peaks indicating 
resonance frequencies. 

The results of this test, when conducted with the soft- 
walled test section and the uniform centerbody are shown in 
Figs. 12 and 13. Figure 12 was obtained in the signal- 
insertion section. Comparison of this result with the 
result shown in Fig. 7 for the hard-walled test section 
shows several interesting features. It is evident that 
peaks corresponding to the hard-walled cut-on frequencies 
are present; however, these response peaks are broader than 
in the hard-walled case. The peaks at 1600-1700 Hz increased 
markedly over that observed with the hard-walled test section 
perhaps as a consequence of reflection at the hard-wall, 
lined-wall junction. There Is also a similar increase in 
the peak centered at 3075 Hz, which is of concern for tests 
that were run at 3150 Hz. A stronger response i's also 
observed, at about 4700 Hz, but this is well above the test 
frequencies and is not of concern for interpretation of the 


test data. 
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Frequency Hz 

Figure 12. Frequency characterization with soft-walled test section 

and uniform centerbody installed; measured in signal insertion 
section. 
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Frequency, Hz 

Frequency characterization of soft-walled test section with uniform 
centerbody installed; measured in lined section. 
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Traces of the system response at downstream locations 
within the test section are shown in Fig- 13- Examination 
of Fig. .13 shows that the pressure amplitude decreases as 
expected as the probe is moved down the soft-walled test 
section farther from the signal-insertion section. This 
. attenuation occurs in all regions, but is most evident in 
the frequency range of 2000-3200 for which the liner was 
designed. Again, this is i-n line with expectations from 
theoretical considerations. 

The response peaks centered at 1600 Hz and 3100 Hz are 
relatively more prominent at the start of -the lined -section 
(Fig.. 13) than they are in the hard signal-insertion section 
(Fig. 12), but they are rapidly attenuated by the liner. 

A similar system' response . curve with the variable-area 
centerbody installed is given in Fig. 14. The results shown 
in Fig. 14 were obtained in the signal-insertion, section. 

The installation of the variable-area centerbody results in 
two. additional resonances in the frequency range 3400-3700 
Hz that correspond .to. cut on of radial modes. There is also 
some shift in the relative strengths of. the standing waves. 
•However, the basic response is quite similar to the response 
observed with the straight centerbody installed. The rela- 
tive response observed at station 1 in the test section with 
the tapered centerbody installed is almost identical to that 
obtained with the straight centerbody installed and is not 
shown. 



3 



Frequency, Hz 

Figure 14. Frequency characterization of soft-walled test section .with variable centerbody installed. 
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These preliminary investigations into standing waves 
indicate that few problems are expected at the frequencies 

840 and 2640 Hz. However, difficulties in comparing the 

? 

test results at 3150 Hz with the theoretical, predictions may 
be encountered since the theoretical model does not account 
for the reflections. 

4.2.3 Determination of Liner Properties 

Determination of the acoustic properties of the liner 
presented some difficulties, even though the effect of the 
backing cavities can be modelled simply by theoretical means 
and the resistance of the facing-sheet material was speci- 
fied as 40 rayls by the supplier. The reactance of the 
facing-sheet material was unknown, and it was readily ap- 
parent that the facing-sheet material is inhomogeneous. In 
addition, the bonding of the facing sheet to the backing 
cavities undoubtably altered the effective resistance of the 
material through blockage of' some pores. Consequently, test 
data from the duct for cases with no mean flow and with a 
straight centerbody installed have been used to determine 
effective properties of the liner section. The liner pro- 
perties established in this manner are. then used in making 
theoretical predictions for variable-area cases both with 
and without mean flow. This procedure is considered accept- 
able since the primary purpose of the tests is to establish 
the ability of the theoretical model to predict the effect 
of axial gradients on the acoustic propagation. 
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The experimental data at station 1 provide the modal 
content of the signal entering the lined section. The 
standard theory for wave propagation in straight ducts 
provides the attenuation and phase changes of each mode for 
an assumed value of the liner impedance. The acoustic 
signal predicted at station 3 in this manner is then com- 
pared with the experimental measurements. The liner imped- 
ance Z is varied to obtain the best fit of the theoretical 
and experimental results at each of the three test frequen- 
cies: 840 Hz, 2650 Hz and 3150 Hz. The admittance 3 = 1/Z 
that is found in this manner for each of the. three test 
frequencies is listed in Table 2 . 

These values of admittance are then related to the 
physical characteristics of the liner to check whether the 
admittances are consistent with the known behavior of 
point-reacting liners. The dimensionless impeda nce Z of a 
point-reacting liner is usually modelled by the formula 

Z = R g (l - iu/uo) + z c (49) 

where R is the dimensionless resistance of the facing 
s 

sheet, wo is the characteristic frequency of the facing 
sheet, and Z is the dimensionless impedance of the backing 
cavities. The impedance of the backing cavities is a known- 
function of the circular frequency w and the cavity depth h; 
the functional form involves the Bessel functions of the 
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Table 2. Liner properties at test frequencies 




.. 

Facing Sheet Properties 

Frequency 

Admittance 

62 

’ R s 

CO 0 

rad/sec 

840 

( 0 . 10 ,- 0 . 26 ) 

1.3 

3 . 0,000 

2640 

( 0 . 77 ,- 0 . 02 ) 

■1.29 

- 25,370 

3150 

( 0 . 6 . 9 , 0 . 12 ) 

l.4l . • 

50,280 


Admittance values calculated using R =1.3 and u 0 = 30,000 
* s 

are 

Frequency 3 2 

840 (. 10 , -. 26 ) 

264.0 (.75, -.095) 

(.75, .13) 


3150 
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first and second kinds of' order zero since the radial varia- 
tion of the pressure in the cavities is given by 1 

p = AJ o (~~) + BYo(^) (50) 

c ^ 

The radial particle velocity v is -then 

v = - "pTc £ AJ S(f^ + BYJ(^)3 (51) 

Requiring the radial velocity to vanish at r = P .2 + h and 

t 

evaluating the dimensionless impedance Z = p/pocv at r = R 2 , 
one obtains the impedance of the backing cavity as 

Z = i J o (~~) Y J [^7 (R 2 + h)] - JJ[^ (R 2 + h)] X 
c c C t* 

Yi(^) / Ji(^)YJ[|(R 2 + h)] 

- JS.C^ (Rz +‘ h)]YJ(^) ( 52 ) 

With Z determined from Eq. (52) , the experimentally-det er- 
c 

mined values of Z were substituted into Eq. (49) and the . 
facing-sheet resistance and characteristic frequency were 
determined. The values of R and w 0 that are found from the 
experimental data are listed in Table 2. The nominal value 
■ of the dimensionless resistance of the facing sheet material 
as supplied by the manufacturer is approximately 1.0. Since 
some reduction of the porosity of the material occurred 
during its bonding to the backing cavities, the resistance 
values of approximately 1.3 at the two lower test frequen- 
cies seem reasonable. Similarly the agreement in the values 
of wo at these two test frequencies is considered satisfac- 
tory. At the test frequency 3150 Hz there is more discre- 



pancy in the predicted liner properties. As noted earlier, 
of the test section and the inability of' the theory to pre- 
dict reflections, cause some doubts as to the adequacy of the 
theory at this test frequency; thus, R = 1.3 and w 0 = 30,000 
rad/sec are taken to be the values of the facing-sheet 
parameters . 

The values of the liner admittance corresponding to R = 
1.3 and w 0 = 30,000 rad/sec have been recalculated and are 
listed in Table 2. The comparison between the experimental 
and calculated pressure amplitudes and phases are shown in 
Figs-. 15, 16, and 17 for each of the three- test frequencies. 
It is readily apparent the circumferential variation of the 
signal could not be modelled- even for this simple case; 
indeed, in the 840 Hz case, the signal has a larger m = 1 
component at station 3 than at station 1. This is entirely 
inconsistent with the standard theory and can be explained 
only in terms of the .observed inhomogeneities of the liner 
material. The soft-walled liner was fabricated in two 
separate halves that were-bqlted together; the results shown 
in these figures strongly indicate that the two halves have 
somewhat different acoustic properties. Thus, the values of 
liner admittance have been selected to provide an average 
fit to the data. Consequently, the comparison between 
theory and experiment for the variable-area cases with and 
without mean flow that are presented in later sections can 
be expected to show only general trends and not the specific 
circumferential variation of the pressure amplitude. 
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Figure 15b. Comparison of the theoretical and experimental pressure 
phases with uniform centerbody installed; M = 0, 

B = .1 - .26i, f = 840Hz. 


Pressure Amplitude, N/M 



Angular Position, degrees 

Fiqure 16a. Comparison of the theoretical and experimental pressure 
amplitudes with uniform centerbody installed; M = 0, 

3 = .75 - . 095i , f = 2640Hz. 
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. ' Fiqure 16b. Comparison of the theoretical and experimental pressure 
phases with uniform centerbody installed; M = 0, 

6 = .75 - . 095 i , f = 2640Hz. 
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Fiqure 17a. Comparison of the theoretical and experimental pressure 
amplitudes with uniform centerbody installed; M = 0, 

3 = ,75 + . 1 3 i , f = 3150Hz. 
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Figure 17b. Comparison of the theoretical and experimental pressure 
phases with uniform centerbody installed; M = 0 
6 = .75 + . 1 3 i , f = 3150Hz. ' 



4.3 Comparison of Experimental and Theoretical Results 

The liner properties estimated in the previous section 
have been used in the theoretical-prediction method described 
in Section 4.1. The acoustic amplitude and phase of each 
mode is determined experimentally at station 1, and the 
theory is used to predict the signal amplitude and phase at 
downstream stations. 

The results for a variable-area centerbody without flow 
are shown in Pigs. 18, 19, and 20 for the frequencies 840 
Hz, 2640 Hz, and 3150 Hz, respectively. The case of 2640 Hz 
yields agreement that is as good as the agreement in the 
straight-walled case; although the circumferential variation 
is not predicted, the overall attenuation rate and the phase 
shift of the 2640 Hz signal are adequately modelled. How- 
ever, the overall attenuation rate is overpredicted for both 
the 840 Hz and 3150 Hz cases, and the phase shift of the 840 
Hz signal is underpredicted. Surprisingly, the results for 
cases with flow that are reported in Section 5 exhibit 
better agreement between theory and experiment. 
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Figure 18a. Comparison of the theoretical and experimental pressure 
amplitudes with the variable centerbody installed; M = 0, 
B • .10 - .261, f = 840Hz. 
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Figure 18b. Comparison of the theoretical and experimental pressure 
phases with the variable centerbody installed; M = 0, 

8 = .10 - . 26i , -f = 840Hz. 
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Figure 19a. Comparison of the theoretical and experimental pressure 

amplitudes with the variable centerbody installed; M' = 0, 
3 = .75 - . 095i , f = 2640 Hz. 
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Figure 19b. Comparison of the theoretical and experimental pressure 
phases with the variable centerbody installed, M = 0, 

3 = .75 - . 095i , f = 2640Hz. 
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Figure 20a. Comparison of the theoretical and experimental pressure 
amplitudes with the variable centerbody installed, M = 0 
3 = .75 + ,13i , f » 3150 Hz. 
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Figure 20b. Comparison of the theoretical and experimental pressure 
phases with the variable centerbody installed; M = 0 
6 = .75 + .1 3i , f = 3TB0 Hz. 



5. Annular Ducts with Plow 


In this section the theory for wave propagation through 
. ducts that carry a mean flow is summarized, the mean-flow 
model used in' the numerical studies is described, and the 
comparison between experimental and theoretical results is 
discussed. The theoretical development closely follows the 
development described in R.ef. 34; thus it is described here 
only in a summary form. The implementation of the analysis 
in a computer code differs in some details from the method 
discussed in Ref.- 3^; however, the changes were introduced 
only to simplify the specification of the mean flow. Since 
the changes have no effect on the accuracy of the calculated 
results, they are not described here, and the reader is 
referred to Ref. 34 for the procedure details. 

5.1 Theoretical Solution 

5.1.1 Problem formulation 

We consider the' transmission and attenuation of acous- 
tic waves in an annular duct haying a slowly-varying cross 
section and carrying a sheared mean flow. The mean flow has 
constant" density ,p, 0 and pressure p 0 . Moreover, the mean 
flow is slightly nonparallel; that is, the mean-velocity 
components along and normal to the duct centerline are 
assumed to have the form U(x l9 r) and eV(x 2 ,r), where x 2 = ex 
is. a slow scale and e is a small dimensionless parameter 
representing, for example, the maximum slope of the wall. 
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We. assume that each flow quantity is the sum of a mean 
part and an acoustic part; that is, we express the flow 
quantities as 

U(xi,r) + EiuCx,r,e,t), eV(xi,r) + eiv(x,r,6,t) , 

(53 ) 

e a w(x,r,0,t), po + e ip(x,r, 0 , t ) , p 0 + £iP(x,r,9,t) 


where £j is a small dimensionless parameter characterizing 
the amplitude of the acoustic wave. Substituting Eq. (53) 
into the Euler equations and subtracting the mean-flow 
quantities, we obtain 


.i£ + o — 

3t 3x 


3p , „ ( 3u 3v 1 9w , Vv 
— + Po( 3l + 3? + F 3? + “ 


- *eV 


3P 

3r 


+ 0(s i ) + O(eei) 


(54)' 


3U i IT ^ 4 . V _ E rv — + u 

3t + U SF + v 3F + F7 3x GLV 8r + U 3 Xl 


,£ - (u I?- + V |§)] + 0 Cei) + 0(££. i) 


Po 


3xi 


3r* 


(55) 


9v. + u BF + i_9P = 

3t 3x . po 3r 


a (W) + 0(e 2 ) + 0(d) 


+ 0 (eei ) 


(56) 


~3w . , TT # 3w , 1 3p 

at ax Por -30 


e(— + v |w ) + 0(e 2 J + o(e a ) 
r or 


+ -0(ee i‘) 


(57) 


■For an i'nviscid, adiabatic ' acoustic disturbance, the equa- 
tion of state gives 

p = pc 2 


(58) 
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where c is the mean speed of sound. The expansion will be 
carried out to 0(e.) in Eqs. (5*0- (57), so that ej must be 
less than or equal to 0{z z ) for consistency. 

The problem formulation is completed by specification 
of boundary conditions: the duct walls are lined with point 
reacting acoustic materials whose specific acoustic admit-, 
tances $j and 32 may vary slow-ly with the axial distance, 

6 1 ■= 3^(xi). For a no-slip mean flow, a requirement of 
continuity of particle displacement gives 

(v .- eR'.u)/p = (-1) J '(3 /.poc) + 0(e 2 ) at v = Rr (59) 

J J J 

where Ri and R 2 are the inner and outer radii of the duct. 

5.1.2 Method of Solution 

Lengths, velocities j time and pressure are made dimen- 
sionless using a characteristic length L, the speed of 
sound e, L/ c and poc 2 as reference values. The method of 
multiple scales determines an asymptotic solution of Eqs. 
(5*0- (58'), subject to. the boundary conditions ( 59 ), of the 
form 

p(x,r,0,t) = p, 0 c 2 F (X'i , r ; e ) exp ( in ) 
v(x,r-, 0 ,t) = cG(x! ,r; e)exp(in) 
u ( x , r , 0:, t ) ■ = cH ( x 1 , , r ; s ) exp (in) 
w(x,r, 0 , 1 ) = c S ( x 1 , r ; e ) exp ( in ) 

where 

3n/3t = - a),- 3q/3x = k 0 (xi),, 3rj/30 = m ( 6 l) 

with constant w. Here w .is the- frequency, m is an integer, 
the real part of k 0 is the wavenumber,- and the imaginary 
part of k 0 is the attenuation rate, a 0 . 
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The modal description of acoustic waves has been employ- 
ed in Eqs. (60). Hence, for each .value of the circumferen- 
tial mode number m there will exist an infinite number of 
•radial modes (including those beyond cutoff) corresponding 
to distinct values of the complex wavenumber ko . For strong 
axial variations of the duct, the radial modes interact as 
they propagate through the duct. However, for the slow 
variations assumed in this investigation the coupling of the 
radial modes can be neglected and the various modes studied 
independently. Thus, Eqs. ( 6 0 ) describe a single arbitrary 
mode. 

The separation of the acoustic disturbance into indepen- 
dent modes also contains an implicit assumption that such 
modes are distinct. In those cases • in which two parallel- 
duct modes have identical eigenvalues, the modes interact to 
produce a linear amplification in the axial direction in 
addition to the attenuation contained in the exp(in) factor 

iig ' • h r 7 

(Tester or Zorumski and Mason ). In a duct whose proper- 
ties vary in the axial direction, the difference between two 
eigenvalues is a slowly-varying function; as the difference 
b.ecomes smaller, the interaction between the two modes, 
which approaches the linear amplification factor as the 
difference vanishes, requires a modification of Eq. (60) to 
account for the additional axial amplitude variations of the 
modes (Section 5.2). However, no theory to account for the 
transition from independent duct modes to the linear ampli- 
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fication interaction is now available, and thus we must 
exclude an analysis of such cases from the current study. 

In' terms of Xi and n, the time and spatial derivatives 

are 

a/at = - w(a/an) 

3/30 = m(3/3n) (62) 

3/3x = k 0 (3/3n) + e(a/3xi) 

The functions F, G, H-, and S are expanded in the form 


p = 

Fo 

(xi ,r) 

+ 

ePi (xi ,r ) 

+ . « . 

G = 

Go 

(xi,r) 

+ 

eGi (x ]5 r) 

+ . . . 

H = 

Ho 

(xi ,r) 

+ 

eHi ( x i ,r) 

+ . . . 

S . = 

So 

(xi ,r ) 

+ 

eSi ( x i , r ) 

+ . * . 


As described in Ref. 34, the expansion of Eqs. (54)- 
(58) according to Eqs'. ( 6 0 ) - ( 6 3 ) leads to the following 
problems for' the pressure functions P 0 and Pi : 

L (F«) = F h % Iv-) + h ( “ 2 - k » - - o («) 

to 2 m 2 

f- 2 - + (-D j+1 ia)0 Po = 0 at r- = R (65) 

J J 

and' - 

L(Pi > = I(P 0j U,-V,k 0 ) ' . (66) 

Mi + (-l)4 +1 i W 0 J F 1 = B ^ ( P o 3 U , V , k o , R j ) at r = R h (67) 
where _ ' 

A 

w = w - kU (68) 

The functionals I and B. contain the effects of the axial 

J 

gradients of the flow; the specific forms of these functions 
may be found in Ref. 34. The slow scale xi appears impli- 
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citly in Eqs. (64) and ( 65 )-. If U, R. , and 3. are inde- 

J J 

pendent of xi, Eqs. (64) and ( 65 ) reduce to those used 
extensive-iy in the literature to analyze the wave propaga- 
•tion in uniform annular ducts carrying sheared mean flow and 
lined with materials having uniform acoustic properties. 

For a given to, U(xi), 3,(xi) and R.(xi), Eqs,. (64) and ( 65 ) 

J J 

can be solved to determine the eigenvalue k 0 (xi) and its 
corresponding eigenfunction ip(r;xi). The solution can be 
•expressed as 

Fo = A(xi )i4(r;xi ) ( 69 ) 

where- A(xi) is still an undetermined function at this level 
of approximation. It is determined by the nature of the 
solution to Eqs. (66) and ( 67 ). Since F 0 and k 0 are solu- 
tions of the homogeneous parts of Eq's. (66) and (67) > the 
inhomogeneous equations (66) and ( 67 ) have a solution if, 
and only if, a solvability condition is satisfied. This' 
solvability condition yields the desired equation for A(xi).' 

To determine this condition, we multiply Eq. (66) by 
rtp(r;xi). Integrating the resulting equation by parts from 
r = Ri to r = R 2 and using the definition of the operator L 
from Eq. (64), we obtain 





3Fi - 
3r 


A 

3r ' 


R] 


Pi 


R] 


R 2 

rli|>dr 

Ri 


(70) 


Using the boundary conditions from Eq. ( 65 ) and ( 67 ) and 
substituting for F 0 from Eq. ( 69 ), one- obtains an equation 
for the amplitude variation in the form 


C~ ~ 3 -' 
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f(xi) + g'( x i ) A = 0 - (71) 

which has the solution 

A(x i ) = A 0 exp[i kidxi] = A 0 exp[ie kidx] (72) 

where 

kr = ig(x i )/f (xi ) (73) 

The functions f(xi) and g(xi) are integrals across the duct 
width of if>,U,V,ko and their derivatives. The integrals 
require numerical evaluation; the specific forms of these 
functions can be found in Ref. 34. Thus the dimensional 
acoustic pressure variation is given by 

p(x,r,e,t.) = A o p o c 2 if>(r;xi)exp i[ (k 0 + eki)dx 

-cot + me] + 0(e) (74) 

where \p(r;xi) and k 0 (xi) are calculated at each station in 
the duct as if the duet walls and mean flow were parallel, 
and 'ki(xi) contains the effects of the axial derivatives of 
the mean flow and duct properties * If the eki term is 
omitted from Eq. (74), the representation of the acousti'C 
signal is referred to as the quasi-parallel approximation 
since it assumes that the acoustic propagation locally is 
the same as in a parallel duct. One objective of this stucly 
is to obtain a comparison of the results of both the -quasi- 
parallel model (i.e., k 0 ) and the ’.non-parallel model, (i.e., 
k 0 + ekj ) with experimental data ,in order to establish the 
validity of the theory. These comparisons are described in 
Section 5-3. 



5.2 Optimum Liner Properties 


Computations based on the method described above have 

attempted to determine the effect of axial gradients on 

"optimum'' liner properties (defined here as those properties 

that produce maximum attenuation of the least attenuated 

il8 q 

mode), since parallel-flow studies } have shown that the 
optimum properties are strongly dependent on parameters such 
as boundary-layer thickness. However, the computed results 
develop a singular behavior as the "optimum 1 ' properties are 
approached. Figure 21 shows results which exhibit such 
behavior; for this specific case, the duct has a constant 
cross-sectional area, the boundary layer is growing, and the 
liner properties are those that a parallel-flow analysis 
identified as "optimum" for the boundary-layer thickness 
that occurs at x = 1.5. Thus the liner is off optimum ex- 
cept at the midpoint where the eigenvalues of two modes 
become equal. The resulting singular behavior has been 
found to be a consequence of a weak nonuniformity in the 
parallel-flow analysis when two modes are equal. The usual 
theory assumes that the modes are independent. However, if 
two eigenvalues k & and k^ are nearly the same, say k^ = k a 
+ A, then 

ik x- ik, x . , , . „ 

p = Ai|> a (r)e + Bi|> b (r)e e (75) 

becomes 
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p = e a (A + Be 1Xl )4(r) + C* (r)e lx ‘ e - lwt+lme (76) 

3, C 

where xj = Ax and i{> is related to the derivative of i{> with 
respect to the eigenvalue. Thus, if the two modes are 
distinct, but close, the acoustic signal has* a slow ampli- 
tude variation with axial distance due to the interaction of 
the two modes in addition to the' slow, non-parallel vari- 
ation included in Eq. (7*0. When the two modes are exactly 

1|6 il7 

equal, Tester and Zorumski and Mason have shown that the 
appropriate acoustic description takes the form 

ik x . 

P = e (A + Bx)l(r) + Cifj (r) ( 77 ) 

cl C 

As "optimum” liner properties are approached and the two 
modes approach one another, the interaction between them 
becomes stronger until it leads to a linear amplification 
factor when the modes are identical. It should be noted that 
this occurs in a very small neighborhood of the optimum ad- 
mittance value. Although the attenuation over a finite length 
section may differ from the "optimal" prediction based on 
a parallel analysis due to the neglect of the linear ampli- 
fication term, the largest attenuation will still be located 
very close to the admittance value calculated without con- 
sidering the mode interaction. This is reflected in the re- 
sults shown in Pig. 21; the local optimum liner properties 
are nearly the same over the entire axial distance 0<x<3, 
yet the singularity develops only very close to the point 
where the liner admittance exactly equals the optimum value. 
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Thus the nonuniformity in the parallel-duct theory is very 

weak; however, in any perturbation analysis, nonuniformities 

become more singular in the higher-order terms — in this 

case in the nonparallel results. The singular behavior 

exhibited in Pig. 21 is seen to approximately the ki = (i/e)/ 

(x-x ) behavior that would be required to obtain the linear 
s 

amplification of Eq. (77) from the slow variation of Eq. $ 
(7^). A result that bridges the singularity in the curve re- 
quires the development of a nonparallel-duct theory that 
accounts for the interaction of nearly-equal modes and thus 
bridges the gap between the classic independent-mode results 
and the linear-amplification results for multiple eigen- 
values. Such an analysis would be analogous to that used by 
50 

Nayfeh to investigate the nonlinear propagation of a wave 
packet in a rigid duct. If we consider only linear, mono- 
chromatic waves, such an analysis would yield an amplitude 
function that is governed by an equation of the form 

h(x) gpr + f(x) + g(x ) A = 0 (78) 

The presence of the d 2 A/dx 2 term, which is generally small, 
becomes important when f(x) approaches sero [f(xi) -*■ 0 is 
responsible for the singular behavior of ki in Eq. (74)] and 
accounts for the linear amplification factor when a multiple 


root occurs. 



5.3 Comparison with Experimental Results 

_ The computer code described in Ref. 34 that implements 

the theory outlined in Section 5.1 has been adapted to the 

flow and wave propagation in the experimental facility 

described in Section 2. The mean flow is assumed to be 

uniform across the duct width except for boundary layers at 

51 

^ the two walls. The integral method of Alber has been used 
to. predict the boundary-layer thicknesses on the two walls 
for the flow conditions at which tests were run. 'The acous- 
tic studies at the three test frequencies of 840 Hz, 2460 Hz 
. and 3150 Hz have shown that the mean boundary-layer thick- 
ness, which is small, has little effect on the acoustic • 
propagation for the test conditions of the experimental 
facility. Thus ■ the refractive effects within the duct are 
primarily due to the axial velocity gradients that occur. 

5.3.1 Straight-wall Cases 

The liner properties used in the theoretical study are 
the same as those estimated in Section 4.2.3 on the basis of 
no-mean- flow results. Comparisons between theory and experi- 
ment have been made for cases in which the duct carries a 
mean flow. Two cases, in which the mean-flow velocity is 
approximately 105 feet per second and 200 feet per second-, 
have been considered for each of the three test frequencies. 

The theoretical results in these cases agree with the' 
experimental data as well as the original no-mean-flow 
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results from which the liner admittance was estimated. 
Several examples are shown in, Figures 22-24. As noted in 
Section 4, the amplitudes of the pressure- disturbance have a 
strong m = 1 component at station 3 that is not predicted by 
the theory, due to the inhomogeneous liner properties. 
However, the -average attenuation of the duct is predicted 
quite well; in the cases that are not shown in Figs. 22-24 
the agreement is as good as for the cases shown, and the 
results for the phase of the acoustic disturbance are ex- 
cellent in all cases. The basic conclusion from these 
results is that the liner admittance is not altered signi- 
ficantly by the presence of the mean flow-. 

5*3-2 Variable-area Cases 

For the variable-area centerbody, the flow in the 
signal-insertion section upstream of the variable-area test 
section was maintained at nominal values of either -105 feet 
per second or 184 feet per second. The flow within the test 
section reached a maximum value of approximately 300 feet 
per second for the faster of the two cases. This keeps the 
Mach number below -0.3 and permits comparison with the incom- 
pressible theoretical model. 

Both the quasi-parallel approximation and the non- 
parallel perturbation approximation were calculated for 
comparison with the experimental data and both theoretical 
curves are shown in those cases in which they differ. 

Figures 25 and 26 show results for 840 Hz at flow speeds of 
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Figure 22. Comparison of the theoretical and experimental pressure amplitudes 

■ with the uniform centerbody install ed; v - 105ft/ sec, 8 = .10 - .26i, 
f = 840Hz. 
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Figure 23. Comparison of the theoretical and experimental pressure amplitudes with 
the uniform centerbody installed; v = 105ft/ sec.,' 3 = .75 - .0951, 
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105 feet per second and 184 feet per second, respectively. 

As before, the liner inhomogeneities generate a strong m = 1 
circumferential-mode component at the third station. However, 
for the lower speed, the average attenuation rate is predicted 
very well by the non-parallel theory in contrast with the 
quasi-parallel theory which overpredicts the attenuation 
rate (Fig. 25a). Both theoretical models give the same 
phase shift (Fig. 25b). At the higher speed, the experi- 
mental pressure amplitude (Fig. 26a) has a somewhat erratic 
circumferential variation that provides little help in vali- 
dating one theoretical result to the exclusion of the other. 
Both theoretical results give good agreement with the mea- 
sured phases, as shown in Fig. 26b. 

At the higher frequencies, the results are also good. 
Figures 27 and 28 present comparisons between theory and 
experiment for flow velocities of 105 feet per second and 
184 feet per second respectively. The pressure amplitudes 
are shown in Figs. 27a and 28a, and the agreement is satis- 
factory in both cases; there is no appreciable difference 
between the two theoretical results. However, the phase 
predictions from the quasi-parallel theory and the non- 
parallel theory are slightly different for 'this case, and in 
both Pigs. 27b and 28b, the non-parallel results are seen to 
be slightly better. 

At the highest frequency of 3150 Hz, the two theore- 
tical results are nearly identical. The comparison with 
experimental data is shown in Figs. 29 and 30. The overall 
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Figure 25b. Comparison of the theoretical and experimental pressure phases with 
the variable centerbody Installed; v = 105ft/sec., f = 840Hz, 

6 ~ .1 - . 26i . 
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Figure 26a. Comparison of the qua si -pa rail el , nonparallel and experimental pressure 
f" 1 - 840Hz $ 1 th6 V 26i ab1S centerbody Installed; v = 184ft/ sec., 
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Figure 26b. Comparison of the theoretical and experimental pressure phases with 
the variable centerbody installed; v '= 184ft/ sec., f = 840Hz, 

A = 1 _ 
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Figure 27a. Comparison of the quasi-parallel, nonparallel and experimental pressure 
2640 Hz* 1 *6 the 7 g aria ^® i centerbod y installed; v = 105ft/sec., 
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Figure 27b. Comparison of the quasi-parallel , nonoarallel and experimental pressure phases 
with the^variable centerbody installed; v = 105ft/sec., f = 2640 Hz, 
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Figure 28a. Comparison of the ouasi-paral lei , nonparallel and experimental pressure 
amplitudes with the variable centerbody installed; v = 184ft/sec., 
f = 2640 Hz, 8 = .75 - .0951. 
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Figure 28b. 


Comparison of the quasi -parall el , nonparallel and experimental 
phases with the variable centerbody installed; v = 184ft/sec 
f = 2640 Hz, 3 = .75 - ,095i. 
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Figure 29a. Comparison of the theoretical and experimental pressure amplitude with 
the variable centerbody installed; v = 105ft/sec., f = 3150 Hz, 

3 - .75 + ,l'3i . 
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Figure 29b. 


Comparison of the theoretical and experimental pressure phase with the 
variable centerbody installed; v = 1 05ft/sec . , f = 3150 Hz. 

6 = .75 + .131 . 
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Figure 30a. 
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• Comparison of the theoretical and experimental pressure amplitude with 
the variable centerbody installed; v = 184ft/sec. s f = 3150Hz, 

p = . /P + i 31 . 
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Fiqure 30b. Comparison of the theoretical and experimental pressure phase with the 
variable centerbody installed; v = 184ft/ sec., f = 3150 Hz, 

3 = .75 + . 1 3 i . 
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attenuation of the amplitudes is again satisfactorily pre 
dieted although the circumferential variation is not. The 
phase results are quite good although there is a fair amount 
of scatter in the experimental results. 

The general agreement between theory and experiment is 
considered satisfactory, discounting the difficulties intro- 
duced by the liner inhomogeneities. These difficulties also 
appeared in the constant— area, no— mean flow cases for which 
the theory is well established and can not be considered as 
a shortcoming of the non-parallel theory. Indeed, the non- 
parallel theory has proven to yield better results than a 
quasi— parallel approximation in those instances in which 
there was a difference between the two approximations. The 
agreement for cases with a me ( an flow that are reported in 

i 

this section has been shown to be better than that^ for cases 
without flow; why this should be so is not clear. 

Overall, the comparisons lend support to- the validity 

of the modeling of the effects of the axial variations that 

occur In the duct. In addition, the need to develop a 

* 

procedure to calculate the effect of circumferential vari- 
ations in the duct properties is clearly present. 



6 . Summary 


An experimental facility has been constructed for the 
study of acoustic wave propagation in annular ducts that 
carry mean flows with both axial and transverse velocity 
gradients. The test section can be selected to have either 
a rigid or lined outer wall; the cross-sectional area can be 
either constant or variable by selection of either a straight 
or variable-radius centerbody, which is rigid. The duct 
liner is point reacting and consists of a Peltmetal facing 
sheet covering cellular backing cavities that are nominally 
tuned to a frequency of approximately 3000 Hz. 

The design of the duct is such' that the mean boundary 
layer is thin within the test section; ambient noise levels 
are controlled by means of two line mufflers and an anechoic 
duct termination so that a good signal to noise ratio is 
achieved. Pour acoustic drivers provide the input signal to 
be studied, and control of the relative amplitudes and 
phases of the drivers permits some control of the complexity 
of the circumferential variation (spinning modal structure) of 
the acoustic disturbance. The duct has been instrumented to 
permit acquisition of narrow-band data for determination of 
the relative amplitudes and phases of the acoustic distur- 
bance along the outer wall of the test section. The test 
sections ahe provided with sufficient microphone probe 
locations to permit resolution of the modal content of the 
signal up to the third circumferential mode. 
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Preliminary experimental tests with both lined and 
hard-walled test sections have established those frequencies 
at which resonances within the duct occur. The resonant 
frequencies, which are not suitable for extensive testing 
for comparison with analytical predictions for propagating 
waves, are shown to be the hard-wall cut-on frequencies; two 
additional resonances occur in the signal-insertion section 
when the lined test section is- installed but apparently do 
not cause problems for study .of propagating waves In the 
lined test section itself.. Three frequencies below cut on 
of the first radial mode were selected for detailed studies 
of. the circumferential and axial variation of the acoustic 
pressure and phase in variable-area ducts and for comparison 
.with analytic predictions. 

Analytical formulas, describing the wave propagation in 
hard-wall and lined annular ducts without a mean flow have 

been developed, and a previously-developed computer code for 

. • * 

prediction of propagation in annular ducts with flow has 
been modified to carry out the computational studies for 
•this investigation. All . analytical methods are based on the 
method of multiple scales for , slowly-varying duct proper- 
ties. Results from the .computer code were used in the 
preliminary design stage of the experimental facility to 

determine a suitable range of values for the area variation 

* * 

of .the duct and. for the liner resistance. 

The computer code has also been used to investigate the 
behavior of the variable-area- effects in the presence of an 
'"optimum" liner • for a variety of flow conditions. It is 
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shown that the non-parallel results develop a singularity at 
a multiple root of the parallel duct eigenequation as a 
consequence of an improper description of the acous.tic- 
propagation process by the parallel-duct theory. 

Comparison of experimental data and theoretical pre- 
dictions in hard-walled duets without mean flow yields 
generally satisfactory results. The phase results at all 
frequencies are good, and a general intensification of the 
acoustic amplitude in the variable-area section is pre- - 
dieted. At the two lower frequencies, 840 Hz and 2640 Hz, 
the circumferential variation of the pressure amplitude is. 
described well by the theoretical predictions, but at the 
highest frequency, 3150 Hz, this variation is inadequately 
predicted . 

For lined ducts without flow, experimental results for 
a constant-area duct are used to establish the liner acous- 
tic properties by comparison with well-establised analytical 
.methods for constant-area ducts. The experimental data show 
an intensification of the m=l circumferential mode in many 
instances' that is though to occur as a consequence of the 
liner being inhomogeneous in its acoustic properties. 
Consequently," the liner properties are determined to be 
those that give the best average fit of the circumferential 
variation, and these liner properties are used to establish 
the ability of the analytic procedures to predict the general 
effects of the area variations for cases both with and • 
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without a mean flow. Tests with flow in a constant-area 
duct have shown that the duct liner properties do not change 
appreciably in the presence of a mean flow. 

For cases without flow, the general agreement between 
theoretical results and experimental data is good at only 
the 2640 Hz frequency. However, the computer code for cases 
with flow predicts variable-area effects that are in rea- 
sonable agreement with the experimental data. Although the 
circumferential variations due to the liner inhomogeneities 
are not predicted, the general effects of the axial varia- 
tions of the mean flow and cross-sectional area are' pre- 
dicted quite well. The non-parallel theory is shown to 
yield better results than a quasi-parallel theory in those 
instances in which they yield different results. Overall 
the comparisons lend support to the validity of the modeling 
of the effects of the axial variations that occur within the 
duct, while pointing to a strong need to develop a procedure 
to determine the effects of circumferential variations in 


the duct properties. 
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